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Abstract 

We investigate relaxation and dephasing of an electron spin confined in a semiconductor quantum dot and subject to spin-orbit 
coupling. Even in vanishing magnetic field, B — 0, slow noise coupling to the electron's orbital degree of freedom leads to dephasing 
of the spin due to a random, in general non-Abelian Berry phase acquired by the spin. For illustration we first present a simple 
quasiclassical description, then consider a model with 2 orbital states only, and finally present a perturbative quantum treatment 
appropriate for an electron in a realistic (roughly parabolic, not too strongly confining) quantum dot. We further compare the effect 
of different sources of noise. While at large magnetic fields phonons dominate the relaxation processes, at low fields electron-hole 
excitations and possibly 1/ f noise may dominate. 
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1. Introduction 

The demonstrations of coherent single-electron spin con- 
trol and measurement [1,2,3] in semiconductor quantum 
dots have opened exciting perspectives for solid state quan- 
tum information processing with spin qubits [4] . More re- 
cent work [5,6,7,8] has revealed further potential of spin 
coherence, which greatly extends the possibilities of next- 
generation spintronic devices. The key behind these emerg- 
ing technologies is the long spin coherence time in semi- 
conductor materials. Spins, unlike orbital electron degrees 
of freedom, do not couple directly to the various sources of 
electric noise present in typical solid-state environments. 

Most of the traditional techniques for addressing and ma- 
nipulating spins in semiconductors have revolved around 
some form of electron spin resonance (ESR), be it through 
external magnetic fields [3] or effective internal ac fields 
based on the spin-orbit interaction. Indeed, spin-orbit in- 
teraction has been proposed theoretically as a way of co- 
herently controlling the spin of confined electrons purely by 
electrical means [9,10,11,12,13], and important experimen- 
tal progress has been made in this direction [14,15]. By the 
same token, it has long been understood [16,17] that spin- 
orbit interaction is one of the main mechanisms by which 
electron spins decay and lose coherence in semiconductor 
heterostructures [18,19,20,21]. 

As we will discuss in this paper, in the particular case of 
an electron confined in a quantum dot, a time-dependent 
(fluctuating or controlled) electric field introduces via the 



spin-orbit coupling a non-Abelian geometric phase (a gen- 
eralization of Berry phases) into the spin evolution. This 
connection between spin-orbit interaction and geometric 
phases has been noted previously in the context of pertur- 
bative analysis of the spin decay of trapped electrons [22]. 
A similar connection had been discussed for free electrons 
in the presence of disorder scattering [23]. 

The geometric character of spin evolution under elec- 
tric fields has striking consequences both for spin-orbit me- 
diated spin relaxation and decoherence as well as for co- 
herent spin manipulation strategies. Geometric spin evo- 
lution under controlled gating is potentially robust, since 
it is not affected by gate timing errors and certain control 
voltage inaccuracies. In the case of spin decay, the non- 
Abelian character of the spin precession under a noisy elec- 
tric environment results in a saturation of spin relaxation 
rates at low magnetic fields [16] through a fourth order (in 
the spin-bath coupling) process previously overlooked in 
the literature [19,24,25]. This spin decay mechanism, which 
can be called geometric dephasing, requires two indepen- 
dent noise sources coupled to two non-commuting compo- 
nents of the electron spin, whereby the non-Abelian prop- 
erties of the SU(2) group become relevant. To second 
order in the spin-bath coupling we also note that a differ- 
ent source of fluctuations other that piezoelectric phonons, 



1 A different phenomenon, also called geometric dephasing, was dis- 
cussed in Ref. [26]. There geometric manipulations of spins in finite 
magnetic fields and the presence of dissipation were considered and 
path-dependent (geometric) contributions to dephasing were found. 



namely electron-hole excitations in the metallic environ- 
ment (ohmic fluctuations), dominate the spin relaxation 
at low magnetic fields. The reason is the higher density of 
ohmic fluctuations at low energies as compared to phonons. 

The present paper is organized as follows. In Sec. 2 we will 
present qualitatively the main concepts and consequences 
of the geometric character of the electrically induced spin 
precession in leading order in the ratio between dot size to 
spin-orbit length, xq/1 80 . In Sec. 3 we consider a model sys- 
tem based on only two orbital states of an electron with 
spin. This model helps to understand the geometrical evo- 
lution of the spin. In Sec. 4 we will perturbatively derive 
the effective Hamiltonian for an electron in a quantum dot 
under a fluctuating electric field taking into account all or- 
bital states. This will allow us to analyze the spin relaxation 
and dephasing under realistic conditions. 

2. Geometric spin precession of a strongly confined 
electron 

Electric fields applied to a quantum dot structure induce 
displacements (and possibly deformations) in the confining 
potential. In the presence of spin-orbit interaction this will 
lead to a peculiar geometric evolution of the spin state of 
the confined electron with important consequences for the 
relaxation and manipulation of the spin. We will approach 
the problem by first considering the spin precession due to a 
geometric phase acquired by the spin of a strongly confined 
electron, when it is adiabatically transported along a given 
trajectory in a 2DEG in the presence of spin-orbit coupling. 

In semiconductor 2D heterostructures the spin-orbit cou- 
pling takes the form (h = 1 throughout this work) 



^SO = OL (fiyCFx - PxO-y) + fl (py&y fix^x) 



(1) 



Here <x/2 and p are the spin and momentum operators, 
while a and (3 are the Rashba and linear Dresselhaus cou- 
plings, which can be lumped into the spin-orbit tensor 




(2) 



It sets the scale for the spin-orbit length l so = ^/\ det A so | = 

[m^/\a 2 — j3 2 \^j . The effective strength of the spin-orbit 

effects in a quantum dot of size xq is in general proportional 
to some power of the ratio xo/l so . In typical GaAs/AlGaAs 
semiconductor heterostructures l so ~ 1 — 5/im, while xq ~ 
30 — lOOnm so that this ratio is usually quite small, of the 
order of 0.02. Other materials, such as In As, have a much 
stronger spin-orbit length, in the l so ~ lOOnm range. 

By classical intuition we can anticipate the main effect. 
We consider an electron in a very strong confinement, forced 
to move along a path C with trajectory Rc(t). Eq. (1) sug- 
gests that the spin-orbit coupling makes the spin precess 




Fig. 1. The geometric precession due to spin orbit interaction for an 
electron adiabatically transported along a path in a 2DEG (left) is 
equivalent to the changing orientation of a sphere rolling on a plane 
along a path (right) which is related by a simple transformation to 
the electron's path. 

under an effective magnetic field B so = ^-pA^ 1 , which 
couples to the spin similar to a Zeeman term except that 
the field depends on the electron's momentum. It raises the 
question as to what 'value' one should use for operator p. 
For a strongly confining potential it turns out that we can 
simply substitute p — > mRc . Hence 



(3) 



From this we derive a spin precession governed by the fol- 
lowing SU(2) operator 



U ad (t) =Texp ^-i jf dtB so • <x 
= Pexp ^— i J dRc A^cr 



(4) 



Here T and P stand for time- and path-ordering operators, 
respectively. The label 'adiabatic' in f7 a d refers to the con- 
straint of slow paths, \Rc\ <C xqujq, typically assumed in 
most works on Berry phases [27]. As is apparent from Eq. 
(4), due the peculiar dependence of B so on the velocity Rc, 
the total "geometric spin precession" for propagation along 
a given path C depends only on the geometry of C itself, 
not on the time dependence of Rc . 

Another line of arguments leading to this result was 
pointed out in Ref. [28]. It is based on the observation that 
H so can be diagonalized to first order in xo/l so by a canon- 
ical transformation exp (—if A" 1 cr), which in turn implies 
that in a small dot the effect of spin-orbit coupling moving 
along a given path can be gauged away by a path-dependent 
gauge transformation /7 a d that rotates the spin just as in 
Eq. (4). 

The evolution operator U a a is a group element in SU(2). 
However, it can also be mapped onto a SO (3) rotation of 
a 3D solid, since both groups are isomorphic up to a sign. 
The natural question arises, what is the 3D rotation cor- 
responding to Z7 a d for a given path? Is there an intuitive 
visualization that tells us how the spin is rotated as the 
containing quantum dot is moved? Remarkably, the SO (3) 
isomorphic form of C/ a d (changing the SU(2) generators fr/2 
by the SO (3) equivalent A) has a very similar form to the 
operator that gives the orientation of a sphere of radius Ro 
that rolls on a plane without slipping or spinning along a 
path C parametrized by Rc(t), 




Fig. 2. The length of a straight path transport of a small dot required 
to perform a pseudospin flip as a function of the angle of the path. 
The three ellipses correspond to different ratios (3 /at but the same 
value of /so- Dotted ellipse: (3 /a = 1/2; dashed ellipse: (3 /a = 1/4, 
solid ellipse: (3 /a = 0. The area of the ellipses remains constant and 
equal to 7r/g . 



U 8ph = Pexp (^-i dR c > A S p h A 




(5) 
(6) 



This picture of the geometric spin-precession in terms of 
a rolling sphere is illustrated in Fig. 1. The radius of the 
sphere is fixed by the spin-orbit length, Rq = l so /2. If only 
Rashba coupling is present there is nothing more to this 
mapping, since in such case A sp h oc A so . However, if Dres- 
selhaus coupling is also present the paths C and C are not 
exactly the same, but are related by a simple 2D rotation 
and scaling transformation Rcf(t)\~^ h = 2Rc(t)X~J . 

This geometrical picture makes it clear that transporting 
the spin along a sufficiently long straight path C in a certain 
direction in the 2DEG will flip the pseudospin. In the case 
of the sphere moving along a straight path C the distance 
needed for a specific rotation does not depend on the direc- 
tion, but for a spin rotation, due to the nontrivial relation 
between C and C when both Rashba and Dresselhaus cou- 
plings are present, the required distance is non-isotropic in 
a 2DEG. In the direction [110] the spin- flip distance can 
be strongly enhanced. This non-isotropic property is repre- 
sented in the Fig. 2 for different ratios of (3/ a. In average, 
the distance that must be covered to induce a pseudospin 
flip is 7r/ so /2, but the variation as a function of the direction 
can be large whenever a approaches (3. 

3. Two orbital states with spin 

3.1. The model 

Before analyzing the dynamics of an electron spin in 
a typical quantum dot we consider a model system with 
merely two orbital levels, |0) and |1), which are occupied 



Fig. 3. Pseudo-magnetic fields h±. 

by a single electron with spin- 1/2. This model is not valid 
to describe electrons confined in (near) parabolic quantum 
dots, where all orbital level are important. Yet, it allows us 
to get a feeling for the geometric phases acquired by the 
pseudo-spin states. 

The Hilbert space of the model system is spanned by 
the four states |0, |), |0, j), |1, |), |1, j). We introduce two 
sets of Pauli matrices: r a for orbital degrees of freedom 
and a a for the spin. Thus, for example, f x \0) = |1), while 
°~x\ T) = | I)- I n the presence of spin-orbit interaction the 
Hamiltonian of the model system reads 



dot 



1 A 

- - er z 
2 



(7) 



where e is the orbital level splitting, while the vector b char- 
acterizes the Rashba and Dresselhaus spin-orbit coupling. 
Choosing real orbital wave functions one can easily check 
that this is the most general form of the spin-orbit coupling 
allowed by time-reversal symmetry. 

Next we introduce the noise which is assumed to couple to 
the orbital degree of freedom. The general form of coupling 
is 



T x X , 



(8) 



where in general Z and X are quantum fields of a bath 
governed by the Hamiltonian i^bath- 

For slow noise, i.e. for noise with power spectrum with 
much weight at low frequencies, instead of considering the 
coupled quantum dynamics of the dot and the bath, it is 
sufficient to treat Z and X as classical stochastic fields, in- 
dependent of each other. Thus, for transparency, we shall 
consider the dynamics of the dot subject to the external 
random fields Z(t) and X(t) and, then, average over real- 
izations of these fields. 

The resulting total Hamiltonian 

11 1 1 

#dot = --ef z - -f y b - & - - Z(t) f z - - X(t) r x . (9) 

shows that the direction of the vector b is the natural quan- 
tization axis for the spin in our problem. Introducing the 
corresponding spin basis states |+)b and | — )& we see that 
the problem factorizes into two subspaces: | \£ orbital) |+) b 




dt 



(e 2 + b 2 fl'- 



Z(t)X(t) 



(13) 



The dots in this relation denote contributions which vanish 
for closed contours and, as can be shown, do not cause 
dephasing. 

Several notes are in order. Since the evolution of the 
pseudospin is described here by a single phase A<l>, rotations 
due to two different loops in the X — Z plane commute. 
Thus we obtain an Abelian Berry phase. Indeed, as we have 
seen, the subspaces + and — do not mix. If, initially, we 
had a superposition a\g+)\+)b + f3\g-)\— )&, the absolute 
values \a\ and \/3\ are conserved and only the relative phase 
changes due to the adiabatic evolution. 

Treating the quantity Z(t)X(t) as a Gaussian stochastic 
field we obtain the dephasing rate as [29] 



Fig. 4. Closed contours traversed by the pseudo-magnetic fields h± . 

and | ^orbital) | — )b- Within each of the two subspaces (±) 
the dynamics is governed by the Hamiltonian 



#4 



1 A 1_ 

-2 6TzT 2 bT ' 



y-\z(t)r z -\x(t)T x 



(10) 



1 



2 (e 2 + b 2 ) 



3 S zx(^ 



0), 



(14) 



where S z j c is the spectral density of ZX. We estimate this 
quantity as 



S zjc (oo = 0) - 2 / duj(j 2 S x (w)Sz(w) , 



(15) 



where b = |6|. We have introduced the pseudo-magnetic 
fields h± acting on the orbital "pseudo-spin" . These fields 
are illustrated in Fig. 3. The fields h± differ only in their 
^/-component, h± iV = ±b. Therefore we obtain \h+\ = \h-\ 
and - consistent with Kramers' theorem - we find that the 
energies of the ground states in the two subspaces, \g+) |+)& 



and \g- 



)b, coincide, E 9i + — E g 



-(l/2)|/i. 



+ ■ 



3.2. Random Berry phase in fluctuating fields 

Next we study what happens when the fields Z(t) and 
X{t) vary in time. In particular, we assume that they tra- 
verse slowly a closed contour in the X — Z plane, as illus- 
trated in Fig. 4. Then the pseudo-spin orbital eigenstates 
states follow their respective pseudo-fields. In addition to 
the dynamical phase each state acquires a Berry phase. As 
one can see from Fig. 4, the Berry phases acquired by the 
two ground states from the subspaces "+" and "-" are op- 
posite in sign 



3>± = ±- / dcp cost 



(ii) 



where the angles (p and are introduced in Figs. 3 and 4. 
Since the dynamical phases for the two ground states are 
the same, the relative phase acquired between the two is 



dip cos t 



(12) 



If we express the angles cp and in terms of the fields 
X and Z and expand assuming X, Z <C e, 6, we find after 
some algebra 



where the limitation of the integration by the temperature 
can be justified by noting that any field can be treated 
classically only at frequencies lower than temperature. A 
fully quantum mechanical analysis confirms this assertion. 

4. Electron spin in a quantum dot 

4.1. Effective Hamiltonian 

Lifting the restriction to two orbital states, we now con- 
sider the full problem of a single electron confined to a lat- 
eral quantum dot by the potential V(r) in the presence of a 
magnetic field B. To be specific, we assume the field to be 
oriented parallel to the plane of the dot, but our procedure 
can be generalized to arbitrary directions [30]. The static 
part of the Hamiltonian then reads 



Ho 



2m* 



V(r) 



9VB 



B-a + H Si 



(16) 



The magnetic field couples to the electron through a Zee- 
man term with material specific q-factorPI. The last term, 
defined in Eq. 1, describes the Dresselhaus ((3) and Rashba 
spin-orbit couplings (a) between the spin cr of the electron 
and its momentum [31]. For a dot of size xo the typical 
energy of the spin-orbit coupling scales as ~ maxa,/3/xo, 
while the level spacing scales as coq ~ l/(ra*#o). Therefore, 



2 For GaAs m* w 0.067 ra e , g « —0.44, longitudinal and transverse 
sound velocities vi = 4.73 • 10 3 m/s, vt = 3.35 • 10 3 ra/s, piezoelectric 



const. h 14 = 1.4 • 10 9 eV/ 



(e*hlJl05(27r)*p(3/vf 



m, 



5.3 • 10 3 Kg/m^ 



density p 

tl4/ iw^/i; u i t ^/ For the quantitative analysis we 
considered Dresselhaus coupling with Xso — l/(jn*(3) = lpm. 



PU(t)P = 
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Fig. 5. Left: Single electron lateral quantum dot in a magnetic field, 
which lifts the ground state degeneracy. Virtual transitions to excited 
states are induced by weak fluctuations of the external fields SE x (t), 
8E y (t). Right: Graphical representation of the evolution operator. 
Virtual transitions to excited states n (wavy lines) are integrated 
out to yield an effective Hamiltonian within the doublet subspace. 

for dots with small level spacing, uoq < IK, the spin-orbit 
coupling cannot be treated perturbatively. 

Next we account for time-dependent fluctuations of the 
electromagnetic field, which add a term 



sv(t)=x»(t) 6, 



(17) 



to the Hamiltonian. (A summation over repeated indices, 
such as /i, is assumed throughout.) The terms M denote 
independent operators in the Hilbert space of the confined 
electron (e.g., x, y, x 2 , . . . ), while the terms X^(t) denote 
the corresponding fluctuating (in general quantum) fields 
(e.g., 5E X , 5E y , V X 5E X1 ...). They may be generated by 
various environments, such as phonons, localized defects, 
or electron-hole excitations. Information about their spe- 
cific properties is contained in the spectral functions, to 
be specified later. Note that this formulation covers also 
quadrupolar fluctuations. 

In case of time-reversal symmetry the ground state of 
the dot is two- fold degenerate. This degeneracy is split in 
an external magnetic field. If q\±bE <C c^o, and as long 
as the noise is adiabatic with respect to the orbital level 
splitting, T <C c^o, the dynamics of the spin remains con- 
strained to these two states. Under these conditions (fol- 
lowing the method described in Ref. [32]) we can derive an 
effective Hamiltonian for the two lowest eigenstates \a = 
±) of Ho by expanding the evolution operator U(t) = 
Texp{— i Jq dt f SVi n t(t')} and projecting to the subspace 
{\cr}} as sketched in Fig. 5. This results in 



PU(t)P = 1 



/' 

Jo 



dhPSViMP 



(18) 



ti>t 2 
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P5V int (t 1 )(l-P)5V int (t 2 )P + ... . 



ti>t 2 



Here 5Vi nt (t) denotes the fluctuating part of the Hamilto- 
nian in the interaction representation and P = \a)(a\. 
We separated terms that involve direct transitions between 
the two lowest states from transitions via excited states. In 



the spirit of an adiabatic approximation, these latter pro- 
cesses can be integrated out to yield an effective Hamil- 
tonian in the two-dimensional subspace. Technically, this 
is performed by introducing slow and fast variables, t = 
(h + t2)/2 and r = t\ — £2, in the last term of Eq. (18), 



-it(e 



,)-ir[i(e„+e a ,)-e»] 6V an (h)6V na , (t 2 ) 



expanding the interaction potential in r as 5V(ti^) ~ 
SV(t) ± J + • • • ? an d integrating with respect to r. 

Here e a and e n denote the eigenenergies of the lowest dou- 
blet and higher eigenstates of i^o, respectively. In this way 
the last term in Eq. (18) becomes local in time. Retaining 
only processes up to 2nd order, we find an effective Hamil- 
tonian within the lowest-energy two-dimensional subspace, 
characterized by the 'pseudospin' Pauli matrices (J x ,y,z, 



if eff = - l -B eS ~a z + X»Cf ■ a + X^C™ ■ a 
+ i (x"X v - X"X v ^j C$ ■ a . 



(19) 



Due to the spin-orbit coupling, which is not assumed to 
be weak, eigenstates do not factorize into orbital and spin 
sectors (hence the term 'pseudospin'). The static effective 
field, £? e ff = (e+ — accounts for the spin-orbit renor- 

malization of the ^-factor and defines the z direction in the 
doublet space. The couplings determining the effec- 

tive fluctuating magnetic fields felt by the pseudospin, are 
given by 



<J,<J f 



a, a' 



M o v , 



(20) 
(21) 
(22) 



The summation is restricted to excited states of higher 
doublets n/^a'. We do not provide explicit expressions 
for the eigenenergies e a , e n , matrix elements 0£ a ,, 0^ n , 
or couplings but below we will evaluate them nu- 

merically and provide quantitative estimates for a generic 
model. We further note that both and C^fJ turn out 
to be transversal to B e ^ , therefore contributing only to re- 
laxation, whereas C^ 2 J has in general also a parallel com- 
ponent that leads to pure dephasing. 

4.2. Geometric phases in B = case 

In time-reversal symmetric situation, (i.e. for 5 = 0), 
the first three terms of Eq. (19) vanish identically [19,30]. 
Only the last term survives, and leads to spin dephasing. It 
has a geometrical origin. To demonstrate this, let us assume 
that the fluctuating (adiabatic) fields are classical. We 
introduce the instantaneous ground states of the Hamilto- 
nian, \<fr n (t)) = |$ n (X /i (t))) defined through the equation 



[Ho + 5V(X^}\$ n (X^) = EniX^niX^)) . 



(23) 



Noting that, to lowest order perturbation theory, the two 
degenerate instantaneous ground states are simply given 
by |<MX M )) * \a) + Zn\ n )( n \SV\a)/(e a - e n ) , we can 
rewrite the last term in Eq. (19) in the familiar form 



(24) 



which shows clearly that the last term is due to a general- 
ized (possibly non-Abelian) Berry phase [27,33,34] acquired 
in a degenerate 2D subspace. In vanishing magnetic field, 
Eq. (24) can be shown to hold to all orders of perturbation 
theory within the adiabatic approximation [30] . If at least 
two linearly independent fluctuating fields couple to the 
dot, they can produce a random Berry phase for the sys- 
tem and cause geometric dephasing at B = 0. When more 
noise components are present, the Berry phase may become 
non-Abelian and all components of the spin may decay. 

4.3. Relaxation and dephasing times 

So far, our treatment has been rather general, applica- 
ble for arbitrary noise properties and dot geometries. In 
its full glory, Eq. (19) describes the motion of the pseudo- 
spin coupled to three fluctuating "magnetic fields" . In gen- 
eral, the dynamics induced by these non-commuting fields 
is complicated. To obtain a qualitative understanding of 
the dynamics we analyze the spin relaxation and pure de- 
phasing times [35], Ti and T 2 * (with 1/T 2 = l/2Ti + 1/T 2 *). 
They are defined only for sufficiently strong effective fields, 
B e fi > 1 /Ti , 1 /T 2 * . In the limit B = we evaluate what we 
call the geometrical dephasing time T geom . For the quanti- 
tative estimate we consider a parabolic confining potential, 



V(r) 



|r 1 2 with level spacing ujq and typical size xq 



1/V^om*. Furthermore, we take into account only dipolar 
fluctuations, X = eSE x xo and Y = eSE y xo coupling to the 
operators Ox = x/xo and Oy = y/xo, respectively. We as- 
sume the two components X and Y to be independent of 
each other, but to have identical noise spectra, Sx(w) = 
Sy(co) = S(u) = 7T£>(a;)cotfi(c<;/2T), with q(uj) being the 
spectral function of the bosonic environment (phonons or 
photons). 

The spectral function q(uj) for phonons can be estimated 
along the lines of Ref. [19]. For the parameters specified in 
Ref. [21] we find for piezoelectric phonons in typical GaAs 
heterostructures at low frequencies, £> p h(^) = #o A p h oo 3 
withAph ~ 4 -10 _6 K _ 2 nm -2 . With these parameters we ob- 
tain relaxation rates generated by the first term in Eq. (19) 
that coincide with those of Ref. [21] at not too high values 
of the field 0. Similar values are obtained for the parame- 
ters of Ref. [36] . For Ohmic fluctuations the spectral func- 
tion is linear at low frequencies, qq(oj) = Xq uj [37,38]. The 



prefactor Xq depends on the dimensionless impedance of 

2 

the circuit, Xq ~ \ Re[Z]. For typical values of the sheet 
resistance of the 2-DEG (10 2 - 10 3 ^/n) we estimate it to 
be in the range 0.1 > Xq > 0.01. For 1/ f noise the power 
spectrum is S(u) = \i/f/\u\. We will further comment on 
its strength below. 

We first estimate the contributions and T 2 ^\ de- 
rived from the three terms (i = 1,2,3) in Eq. (19), for anon- 
vanishing in-plane magnetic field, £? e ff ^> 1/Ti, 1/T 2 *. The 
coupling turns out to be perpendicular to £? e ff [21], 
and for low magnetic fields and weak spin-orbit coupling is 
proportional to |C^| ~ max{a,/3}. This fluctuating 
field therefore contributes to the Ti energy relaxation only, 



(i 



2 |C 



(1)|2 
X I 



t(l)|2 



Sx(B eS ) . 



(25) 



It scales as 1/T^ ~ B 2 max{£?,T} for Ohmic dissipation 
and as ~ B 4 max{£?, T} for phonons. As a consequence, for 
dots with level spacing in the range ujq « 1 . . . 10 K Ohmic 
fluctuations dominate over phonons for low fields with B < 
1...3T. 



The second term, 



max{a 2 ,/? 2 } gives rise 



to both relaxation and dephasing. The two rates 



are 



n 



.(2,-L) .2 
'XX, s\ 



^(2,-L)i2 
y YY,s I 



*(2) 



= 4(|C^J 2 - 
Sxy(u) = | J 



1^(2,11) ,2 
\^YY,s\ 



2|C^| 2 )^y(5 eff ), 
£| 2 )^y(0), 



(2 

XY, 



duo 



1 - cosh(cD/2T)/cosh(cc;/2T) 



with Cj^y denoting the symmetrized component of C^ 2 J 



perpendicular /parallel to T? e ff • Thus, for Ohmic dissipation 
1/T^ vanishes as ~ B 2 max{£? 3 , T 3 }, while for phonons 
it scales as ~ B 2 max{£? 7 , T 7 }. 

is also perpendicular to T? e ff • Its contribution to the 
relaxation is 



1 



(3) 



: 2 l^xy,al 2 ^xy-xy(^eff) , 



(26) 



J XY-XY 



duj 



1 - cosh(cD/2T)/cosh(cj/2T) ' 



with C^} a being the anti-symmetrized component of 
Cffl. Most importantly, with \C^y\ ~ ^max{a 2 ,/3 2 }, 

the rate 1/T-j^ approaches a non-zero value at low fields, 
1/Ti,l/T 2 * <C B <C T, and scales as - max{£ 5 ,T 5 } for 
Ohmic dissipation and ~ max{£? 9 , T 9 } for phonons. 

Finally, at B = the geometric dephasing rate is given 
by l/T-p\ extrapolated to zero field 



3 For larger frequencies the approximation Q p h(u) = Xq uj 3 is not 
valid since the wavelength of relevant phonons becomes comparable 
to the size of the dot. 



4 Note that there is no discrepancy with the scaling max{a, (3} 
quoted in Ref. [19], since they include higher (multipolar) spin- 
flipping phonon contributions neglected here. 
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Fig. 6. Spin relaxation rates for a GaAs quantum dot with level spac- 
ing ooo = IK as function of the Zeeman field. We chose a tempera- 
ture T = 100 mK and Ohmic coupling X n = 0.05. Below B* w 1 T 
spin relaxation is dominated by coupling to Ohmic fluctuations. 
For B < B** ^ 15 mT geometrical spin relaxation due to cou- 
pling to Ohmic fluctuations dominates. For all B values plotted, the 
Bloch-Redfield consistency requirement, B e ff ^> is satisfied. 

= 2\cf Y J 2 Sx Y - X y{u ~ 0) . (27) 

geom 

In our example with only two noise components this process 
dephases only the components of the spin perpendicular to 

The relaxation rates corresponding to the different terms 
in Eq. (19) and various noise sources are shown in Fig. 6. 
Clearly, for external fields 5<5*^1...3T Ohmic fluc- 
tuations provide the leading relaxation mechanism. The 
crossover field B* is not very sensitive to the specific value 
of the spectral parameter Xq and is independent of the 
spin-orbit coupling. Below a second crossover field, w 
15mT, the geometric dephasing induced by Ohmic fluctu- 
ations starts to dominate. This second crossover scale is 
very sensitive to the spin-orbit coupling and temperature, 
scaling as 5** ~ max{a, xq)(T /ujq) 2 . E.g. for a level 
spacing c^o ~ IK and temperature T = 100 mK (T = 
50 mK) the Berry phase mechanism gives a relaxation time 
of the order of 700 fis (20 ms). For even lower temperatures 
or smaller dots with level spacing loq ~ 10K the B — > re- 
laxation time is quickly pushed up to the range of seconds. 

Finally, we comment on the effect of 1/f noise. In 
most cases, the non-symmetrized correlators for 1/f 
noise, needed to calculate correlators as Sxy or Sxx, 
are not known. Yet, for <C T we can provide an 

T 

estimate Sxy(w) ~ J ^Sx(u — £o)Sy(&), and simi- 

-T 

larly for S xy _ X y- The B = geometrical dephasing 
rate due to the 1/f noise can be estimated as T~ e l m w 
\^xy \ 2 ^i/f(T) lo c , where u c is the upper frequency cut-off 
for the 1/f noise [39]. Accounting for the high-frequency 
(Ohmic) noise, sometimes observed to be associated with 
the 1/f noise [40,41,42,43], the estimate becomes T^ m w 

\ C XY\ 2x l/f( T ) T - While the 1 /f noise of background 
charge fluctuations is well studied in mesoscopic systems, 



the amplitude of the 1/f noise of the electric field in quan- 
tum dot systems is yet to be determined. If we assume that 
this noise is due to two-level systems at the interfaces of 
the top gate electrodes, we conclude that in the parameter 
range explored here, the effect of 1/f noise is less impor- 
tant than that of Ohmic fluctuations. However, in quantum 
dots with large level spacings in the low-temperature and 
low-field regime, these fluctuations could dominate over 
the effect of Ohmic fluctuations and eventually determine 
the spin relaxation time. 

5. Conclusion 

We have shown how spin-orbit interaction in single elec- 
tron quantum dots induce pseudospin precession within 
Kramers' doublets when the quantum dot is adiabatically 
shifted along a path by electric fields. The precession de- 
pends solely on the geometry of the path. We have analyzed 
the resulting spin decoherence within a model restricted 
to two orbital states and showed how random geometric 
phases appear in this model. Then we have studied spin 
relaxation and dephasing in a quantum dot considering all 
orbital states. We derived an effective Hamiltonian for the 
lowest Kramers' doublet. The geometrical spin precession 
induced by at least two independent electric fields occurs 
around non-commuting axis for different paths, so that it 
becomes a strictly non-Abelian evolution. This has marked 
consequences for the spin decoherence due to electric field 
fluctuations. In particular the spin decay rate remains fi- 
nite as the external magnetic field vanishes. We estimated 
the rates of relaxation processes due to the coupling to 
phonons and to the ohmic environment. We find that, ow- 
ing to the higher excitation density of ohmic fluctuations 
at low frequencies as compared to phonon fluctuations, the 
former dominate over the latter in the low field relaxation 
rates and provide most likely the leading relaxation chan- 
nel. The underlying spin precession mechanism could also 
be exploited to electrically manipulate spins, namely if by 
applying electric fields the electrons can be transported 
around ~ /im-sized loops. 
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